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A SIMPLE NUMERICAL METHOD J’OR jTfiE -CALCULATION OF 
THE LAMINAR BOUNDARY LAYER 
By K. Schroder 

April 1952 


Every effort is made by the NACA to insure that its published 
translations are accurate reproductions of the original work of 
the authors. Papers are selected for translation on the basis of 
interest and probable usefulness and, although an examination is 
made for technical soundness, the Committee cannot assume respon- 
sibility for the accuracy of the detailed results presented by 
the author in the original paper. The Committee will, of course, 
call attention to any errors observed at the time of the publica- 
tion or subsequent thereto. 

A recent application of equation 21 on page l4 of TM 1317 
indicated that the denominator of the last term, which is given 
as 2 in the NACA translation and in the original German document , 
should be l . 
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A SIMPLE NUMERICAL METHOD FOR THE CALCULATION OF 
THE LAMINAR BOUNDARY LAYER 
By K. Schroder 

April 1952 


In the errata no. 1 issued on this paper, the last term of equa- 
tion (2l) was corrected as follows: 

U(6p+l) - U^p-p) 

Subsequent consideration of this errata indicates that this corrected 

1 2 

term is valid only for the case when k = Therefore, when equa- 

tion (2l) is used for general applications, the last term therein 
should he 
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TECHNICAL MEMORANDUM 1317 


A SIMPLE NUMERICAL METHOD FOR THE CALCULATION 
OF THE LAMINAR BOUNDARY LAYER* 

By K. Schroder 


ABSTRACT 


A method is described which permits an arbitrarily accurate calcu- 
lation of the laminar boundary layer with the aid of a difference cal- 
culation. The advantage of this method is twofold. Starting from 
Prandtl's boundary-layer equation and the natural boundary conditions, 
nothing needs to be neglected or assumed, and not too much time is 
required for the calculation of a boundary- layer profile development. 

So far, the method has been tested successfully in the continuation of 
the Blasius profile on the flat plate, on the circular cylinder inves- 
tigated by Hiemenz, and on an elliptical cylinder of fineness ratio 1:4. 
Above all, this method offers for the first time a possibility of con- 
trol by comparison of methods known so far, all of which are burdened 
with more or less decisive presuppositions. 


OUTLINE 


I. INTRODUCTION 

II. GENERAL REPRESENTATION OF THE METHOD 

III. PRACTICAL EXECUTION 

IV. NUMERICAL EXAMPLES AND RESULTS 

V. REMARKS REGARDING THE CONVERGENCE OF THE ITERATION PROCESS 

*"Ein einf aches numerisches Verfahren zur Berechnung der laminaren 
Grenzschicht." Zentrale fur wissenschaftliches Berichtswesen der Luft- 
fahrtforschung des Generalluftzeugmeisters (ZWB ) Forschungsbericht 
Nr. 1741, Ber lin-Adler shof , February 25, 1943. 
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I. INTRODUCTION 


The flow processes in the laminar "boundary layer may he described 
by Prandtl's boundary-layer equation. If one limits oneself to the two- 
dimensional steady case and introduces, in a suitable region around a 
profile contour C situated in a flow, a curvilinear coordinate system 
s,n, the coordinate lines of which consist of parallel curves and nor- 
mals of C, that equation reads 


3v<, 

v s —2. + v n 
os 



p' (s) 


= 1 ° 2y s 
R 3n 2 


( 1 ) 


when v s , v n signify the velocity components in the s,n system, R 
the Reynolds number, and p = p(s) the pressure distribution along C 
taken from a measurement or calculation 1 . Equation (l) is complemented 
by the continuity equation 


3v s 

bs 


+ 


Svn 

3n 


= 0 


The transformation 

T) = n l/5 Vti = v n J/r 

yields, instead of equations (l) and (2), the equation system 


( 2 ) 


v s 


qy s 

3s 


+ V- 


- a ^ + P'(B) 

OR 


3 2 v s 

W 


(3) 


+ 'tl = 0 

3b 3rj 

in which R no longer appears explicitly. 


w 


^A mathematically complete derivation of equations (l) and ( 2 ) 
based on physically plausible assumptions may be found in H. Schmidt's 
and K. Schroder's report entitled "Die Prandtlsche Grenzschichtgleichung 
als asymptotische Naherung der Navier-Stokesschen Differentialgleichungen 
bei unbegrenzt wachsender Reynoldsscher Kennzahl" (Prandtl's boundary- 
layer equation as an asymptotic approximation of Navier -Stokes' differ- 
ential equations for indefinitely increasing Reynolds number) Deutsche 
Mathematik, 6, Heft 4/5, pp. 307-322. A survey of related literature 
is given in H. Schmidt's and K. Schro’der's report "Laminare Grenz schichten 
I. Teil" (Laminar boundary layers, part I) Luftfahrtforschung 19, 

Lieferung 3, 1942. 
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The following "boundary conditions for the integration of equa- 
tions ( 3 ) and (^) are usually selected in houndary-layer theory as the 
natural ones from the physical point of view. For an initial value s = s Q 
an entrance profile 

v s = v s (s 0 ,T}) 

is prescribed as a function of T} (entrance condition). Furthermore, 
in consequence of the adherence 'of the fluid to the contour, the rela- 
tions 

Ho ’ 0 H 0 ’° 

which are to be interpreted as limiting processes, are to be valid along 
C (adherence condition). Finally, for s-values larger than or equal 

to s 0 the velocity component v s is to converge for 3 > 00 toward 

the velocity U(s) which is connected with the prescribed pressure 
distribution p(s) by 


U( 3 ) U' ( s) = - p’ (s) 


( 5 ) 


(transitional condition). 

The general significance of these boundary conditions will be dis- 
cussed more thoroughly in the second part of the Luftfahrtforschung 
report quoted in footnote 1. Here we shall only point out that the tran- 
sitional condition formulated for 3 — > °° must not be confused with a 
condition for n > °° since for the latter limiting process the veloc- 

ity components converge toward those of the basic flow^. The limiting 

process rj * denotes, on the contrary, the asymptotic transition 

to the boundary values, resulting along C of the outer potential flow 
obtained for R ' — > «> . This can best be made clear by the example of 
the stagnation-point flow at the flat plate, treated in the second report 
indicated in footnote 1 (by the author and H. Schmidt).. Whereas the 
quantity 5 there specified as boundary- layer thickness tends 

like i/l/z toward zero, a quantity d tending toward zero, for 

instance, like 1/ (jS", can be prescribed in such a manner that the flow 
outside of a layer of the thickness d adhering to the contour for 
R — > 00 converges toward the outer potential flow. However, to the 
asymptotic transition toward the boundary values of this assumed potential 


%t is assumed, of course, that this limiting process is meaning- 


ful. 



4 


NACA TM 1317 


flow along C then there corresponds the limiting process 

lim = lim d\/R = <» 

R — »°° R — *°° 

So far, an appropriate existence and uniqueness theorem for this 
boundary- value problem does not exist. However, the results obtained 
with the new method described below show that the statement of the 
problem is perfectly sensible. 

In the literature it has been pointed out more than once3 that for- 
mal power series developments of the function representing the solution 
with respect to t| make the fact plausible that the entrance profile 
cannot be selected completely arbitrarily, but that it is dependent on 
the pressure distribution p = p(s). 

Our method for the determination of the velocity profiles yields a 
numerical solution of the mentioned boundary- value problem with the aid 
of the difference calculation; it is superior to other methods because 
it requires no assumptions beyond equations ( 3 ) and (4) and the boundary 
conditions. In our method, the boundary-layer bonds of the entrance 
profile do not appear directly and thus do not cause any difficulties in 
the numerical calculation. A severe violation of these bonds causes, 
in our method, the variation of the successive boundary-layer profiles 
to become completely disordered. Small violations of these. bonds, in 
contrast, do not exert any considerable effect on the further develop- 
ment of the profile^. 


^Compare S. Goldstein "Concerning Some Solutions of the Boundary- 
Layer Equations in Hydrodynamics," Proc. Cambridge Phil. Soc. 2 6 , 1930, 
pp. 1-30, L. Prandtl, "Zur Berechnung der Grenzschichten" (Concerning 
Calculation of the Boundary Layers) ZAMM. 18, 1938, pp. 77-82, (NACA 
TM 959) and H. Gortler, "Weiterentwicklung eines Grenzschichtprofils 
bei gegebenem Druckverlauf" (Further development of a boundary-layer 
profile for prescribed pressure' variation) ZAMM. 19, 1939, PP. 129-140. 

**L. Prandtl and H. Gortler (reports quoted in footnote 3) 
at the same conclusion, although on another basis. 


arrive 
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II. GENERAL REPRESENTATION OF THE METHOD 


If one introduces into, equation (3) ., instead of s, the new inde- 
pendent variable 


Js 0 U(t) 

under the assumption that U(s) ^ 0 for s = s Q whereby 

dl _ 1 

<3.s U(s) 

is valid, and if one uses the new designations 

u(| ,Tj) = v s (s,t)), U(!) = U(s(|) ), U*(I,T]) = u(! ,T|) - U(l ) 
there follows from equations ( 3 ) and (4) by way of 


Vn AH drj 

1 0(5) J 0 85 

with equation ( 5 ) taken into consideration, our initial equation 


Su* _ S 2 u* _ 1 Su 

SI St) 2 u(|) Sq 


■* Su 


UP _ —ii — 

SI u( 1) s^ 


According to the statement of the problem in the introduction, we 
have to find a solution 


U* = U* ( |,Tj ) 

of equation (7) for all points (|,ri) in the right upper quadrant of 
the plane of the rectangular Cartesian |,rj coordinates? which in 
approaching the straight lines 

|=0 or r| = 0 respectively 


5lf separation phenomena appear, the solution will, in general, be 
of interest only up to the separation point or possibly a little way 
beyond it. 
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tends toward prescribed functions: 


lim = v s (s 0 ,T]) - U( s o) (n - °) (®) 

I — *0 


or 


lim u*(|,T)) = - U(0 (I = 0) (9) 

T|— ^0 

and which vanishes for q — ><» 


lim u*(S,q) =0 (6^0) (10) 

T) 

The fundamental formulation of our method consists in using the 
functional relation (7) - in the sense of the known method of successive 
approximations - for the calculation from a prescribed approximate solu- 
tion which already satisfies the indicated boundary conditions of a 
sequence of corrected functions .which converges toward the actual solu- 
tion of the problem; one substitutes the last obtained approximate solu- 
tion every time on the left side of equation (j) and integrates the 
resulting partial differential equation of the type of the inhomogeneous 
heat conduction equation. 

The examples so far calculated numerically showed that the iteration 
process is obviously convergent. Nevertheless, a general proof of this 
fact would be very desirable and we reserve returning, in a given case, 
to a mathematical examination of . these problems. (Compare also Section V.) 

One may characterize the method by stating as the desired result a 
continual improvement of a given approximate solution in the sense of 
Oseen's method of linearization. Then this linearization of the hydro- 
dynamic equations of motion (which, of course, for the boundary-layer 
flow taken by itself is not permissible) consists in introducing the 
velocity loss u* and in neglecting all nonlinear terms in u*,v and 
their derivatives. From equation ( 3 ) one would thereby obtain 


U(s) 


du* _ cku + u *dU _ 0 

& 3q 2 48 


thus on the left side (aside from the term u*4^ which, however, does 

as 

not alter the character of the equation) precisely the expression which 
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also appears on the left side of our Initial equation (7). 

In integrating (under the houndary conditions (8), ( 9 ), and (10)) 
the differential equation 


chi* _ S2u* _ 
81 8t) 2 


= f(S,n) 


( 11 ) 


into which had heen introduced for abbreviation the function 


"■n 




dn - -J&. 

U(|) ^ 


( 12 ) 


to be regarded as known in the sense of our approximations, one may now 
use successfully the difference calculation. For the homogeneous equa- 
tion this has been done, simultaneously with a proof of convergence, by 
R. Courant, K. Friedrichs, and H. Lewy°. For the inhomogeneous equation 
here dealt with, the proof of convergence together with a formula for 
error estimation may be found in a paper by L. CollatzT. 


If one covers the right upper quadrant of the |,tj plane by a net 
of lattice points with the coordinates 


Ip = pk ha = crZ (p,a = Oj integers) 

(compare fig. l) and introduces at the same time, with a view to later 
applications, the new designations 


u p >0 


u(£p,ha)> u *p,a = u*U p ,T) a ) 


8u 


8u 

C 

-1 


_8I_ 


81 

PyO [2 




^R. Courant, K. Friedrichs, and H. Lewy: "Uber die partiellen 

Differenzengleichungen der mathemat i s chen Physik" (On the partial differ- 
ence equations of mathematical physics). Math. Ann. 100, 1928, pp. 32-7^, 
particularly pp. ^ 7 - 52 . 

Tl. Collatz: "Das Differenzenverfahren mit hoherer Approximation 

fur lineare Differenzengleichungen" (The difference method with higher 
approximation for linear difference equations), Schriften des Math. Sem. 
u.d. Inst. f. agewandte Math.d.Univ. Berlin, Bd. 3 } Heft 1, 1935* 
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there corresponds to the differential equation (11) the difference equa- 
tion of first approximation 


u*p+l,cr - u*p,cr u *p,cr+l _ 2u* p ^ a + 


= f 


p,a 


(13) 


If one selects the step magnitudes k and 2 in I and q direction 
not independent of each other but so that 



( 1 ^) 


equation ( 13 ) is transformed into the simpler difference equation 


u 


p+l,CT = 


u *p,cr+l + u *p,o--l 


+ kf 


p,a 


(15) 


It can be shown that the solution of equation (15) for the corresponding 
boundary-value problem for 2 — >0 and therewith also for k — *-0 con- 
verges toward the known solution of the boundary value problem of equa- 
tion (11). 

Since the values 


u *p,0 (P = 0) u *o,cr ( a = 0) f pf0 (p ^ 0 ,a > 0) 

are known, one may, according to equation ( 15 ) , successively calculate 
all values 



(p ^ o,a ^0) 


progressing stepwise from lattice point column to lattice point col umn . 

Actually, however, we apply another correction at every step in order 
to compensate the systematic error originating by the fact that the 
derivative appearing on the left side of equation (11) 
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was replaced by the difference quotient of first approximation 

u *p+l,cr - u * p f a 
"k~ 

(Compare the following section. ) 

One notes that due to the transitional condition (10) for the entrance 
profile u*q Q , necessarily must vanish for a — and that fp a 

likewise vanishes for a — ^ OT , since even the approximate solution used 
for the formation of f p a was supposed to satisfy the condition (10); 

hence one recognizes that the corrected solution (obtained with the aid 
of the difference calculation in the manner described above) also satis- 
fies' the transitional condition (10). 


III. PRACTICAL EXECUTION 


In practice one may vary the method in such a manner that one does 
not at all require an approximate solution prescribed at the outset in 
the first quadrant of the £,T] -plane; one rather determines this approx- 
imate solution for every step and then improves it to the desired accuracy 
before passing on to the next step. Thus one applies a combined system 
of continuation and correction. 

If one deals with the flow about a profile contour, the initial pro- 
file at the point s = s Q is best taken from the well-known power-series 
developments by Blasius-HIemenz, the coefficients of which for the first 
three terms were given in table form by Howarth^. For reasons of con- 
vergence, these broken-off series will represent a good approximation of 
the solution of the boundary-layer equation only at a small distance from 
the forward stagnation point (s = 0) of the outer potential flow. In 
the permissible range they represent, as it were, an improved stagnation 
point flow. 

Our calculations so far have shown that the series are serviceable 
up to s-values for which the "first boundary layer bond" 

^1 =p'( 3 )=-MS1 (16) 


A 

°Compare L. Howarth: On the calculation of steady flow in the 

boundary layer near the surface of a cylinder in a stream. R & M no. 

1632, 1934. 



10 


NACA TM 1317 


which is a direct result of equations ( 3 ) and (6) is satisfied with 
sufficient accuracy. 

In practice, one has therefore to start the calculation by approxi- 
mating the function U(s) for small s-values as well as possible by a 
polynomial of the form 


U(s) = uqs + u^s^ + UtjS^ 

for the case of a profile symmetrical in free-stream direction, or 
respectively, of the form 


U(s) = u^s + U 2 S ^ + U 3 S ^ 

for the case of a profile unsymmetrical in free-stream direction; one 
may sometimes get by with only two terms. 

After having determined, in the manner described above, the value 
s Q > 0 at which the continuation method may start, one first sets up 
the connection (given by equation (6)) 



(s > s Q ) 


by evaluating the integral on the right side, for instance according to 
the trapezoidal rule. One graphically represents the functions \ = £(s) 
and U = U(s) in a common diagram so that U = U(|) can be immediately 
taken from it. 


The step magnitudes k and Z, connected by equation (l4), must be 
selected so that, first, a sufficient number of subdivision points are 
distributed over the profiles to be calculated, and second, a sufficiently 
rapid continuation in | direction is possible. When profiles of not 
too pronounced S-shape (near the separation point) are to be calculated, 
eight to ten equidistant subdivision points generally will be sufficient 
to define the profile. In upward direction (that is, for large t] values) 
one will have to take so many subdivision points that the profile dies 
out sufficiently gradually toward the asymptotic value U. This provides 
a first indication for the selection of Z and therewith also of k. 

It should finally be remarked regarding the step magnitude k that it 
must be at least large enough to make, for fixed . I and variable tj, 

the derivatives ^ (obtained in first approximation by formation of 
o| 


difference quotients) take a reasonably regular course (compare the 
following discussion) . Hence the lower limit is set for k and there- 
with also for Z . 
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On the other hand, one will he forced to choose the smallest possible 
step magnitude k at points I where the curves u = u(£, const) exhibit 
great curvatures (which occurs particularly directly ahead of the separa- 
tion point), in order to make a sufficiently exact calculation of the 
profiles possible. There Z, too, will necessarily be small. Since, 
however, the boundary- layer thickness has greatly increased at the sepa- 
ration point, one will have there a great many subdivision points dis- 
tributed over the profile. This is in one respect convenient - the posi- 
tion of the separation point is better defined. On the other hand, the 
expenditure of work increases at such points. However, at the end of 
this section we shall point out a possibility of reducing the steps 
in I direction without necessarily having to accept a step reduction 
in T] direction. At the same time we shall then be able to indicate a 
criterion by which the necessity of a step reduction in i direction 
may be recognized. 


Once a certain selection of step magnitudes has been decided upon, 
it is a question of obtaining a first approximation for the values fQ a 
appearing in equation ( 15 ) , in order to be able to execute the first step 
in £ direction. It should be noted that together with the initial pro- 
file at s = s Q also the values of u for values s < s Q may be. taken 
from the series developments. Particularly the values u 1 ^ (that is, 

the profile one step ahead of the initial profile) are thus known. 


We then put for a first approximation of the 


w 


in f. 


0 ,cr- 


d£ 


0 , cr 


oc curing 


c)u 

SI 


= T* 0 ’* " U ~ 1 ' g 

0,cr k 


Therewith 

experience has shown that this Integration may be very conveniently carried 
out with sufficient accuracy by use of the trapezoidal rule with the aid 
of the present subdivision; this can be done purely schematically by 
calculation according to tables. For at the | points where the deriva- 


p 

du 

'0 

as_ 


dr] too can be evaluated numerically. Our calculation 


tives 


cju 

SI 


become very large - whereby the values 



come to 


be of great importance in the calculation of the profiles and must be 
determined relatively exactly as for Instance in the neighborhood of the 
separation point - it will be necessary to select small k (and there- 
with also Z) values so that a sufficient number of subdivision points 
are distributed over the profile to allow application of the trapezoidal 
rule with sufficient accuracy. 
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If one puts, furthermore, with good approximation 


&11 


JO ,0 


U 0 , g +1 ~ ^, 0-1 

2 1 


( 17 ) 


u 1 a and u-^ a may be calculated in first approximation. 


* * 1 1 
u ljl,a and ^1)1,0- 


The values thus obtained will be denoted by 

It will be best to- arrange the entire calculation procedure in the form, 
of a table (compare table I on page 3 l)- With the values obtained 

one will form corrected values of the derivatives 

II, a 

according to the scheme 


H: 


Si 


0,a 


- 

OU]_ 

■ A 

ll ko - U -l,a 

bl 

0 

q 

1 

1 

2 k 


whereupon one obtains (with the aid of table II on page 31) a second 
approximation jugjj a ^ or va l ues u ]_ a with the values 7 a, e a, 


and C a taken from the first table. Whereas the derivatives 


Si 


0 ,CT 


formed in first approximation might show at a few points a an irregular 
course, this will generally no longer be the case for the corrected deriva- 

rsui 


tives 


Su 

S^Jo,a 


The columns for the quantities 






and 


^2 1 a occurring further on in table II will be explained only later. 


This procedure is continued until the values obtained in the third- 
from-last column of the table no longer vary in the desired decimal. In 
the examples we calculated the iteration was carried so far that for every 
step | a the values Up a no longer varied except for an error of 

about l/b to l/2 percent of the maximum velocity U(|p) in each case. 

For the selected step magnitude k this was the case after two to three 
iterations. 


Due to the favorable position of the errors, the profiles calculated 
in the manner described generally show a very smooth course. If the 
u = u(£, const) are concave in respect to the I axis, as is the case 
for instance in the flow about the circular cylinder or the ellipse near 
the separation point (compare fig. 7 and fig. 12), the convergence occurs 
only on one side in the direction from larger to smaller values for u. 
The opposite behavior exists when the course of this curve is convex with 
respect to the % axis as is the case for instance in the boundary- layer 
flow at the flat plate. 
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If one wants to obtain with the described procedure a calculation 
of the u variation as accurate as possible without selecting too small 
a step magnitude k, thereby increasing too much the expenditure in cal- 
culation^ one will find it necessary (as mentioned before) to make at 
every step a correction which takes the fact into account that in setting 
up the basic equation ( 15 ) the difference quotient of first approximation 


only was substituted for the derivative 


ki 

ai 


p,cr 


If one were to select instead the representation of higher approxi- 
mation 




u p+1, cr 


u p-l,cr 


2k 


one would obtain by maintaining equation (14) 


* * * * 
u p+1 , 0 - u p— 1, o' + u p,cn-l + u p,o-l 



+ 2kf 


p,o 


(18) 


instead of equation (15). 

Since this relation, however, (as can be seen immediately) behaves 
considerably less favorably regarding propagation of errors than equa- 
tion ( 15 ), the profiles calculated with its aid will no longer show the 
smooth course mentioned before. Calculation practice has shown that one 
obtains very smooth curves if one writes instead of equation ( 18 ) 


* * 01 

u „ = u_t „ + 2k 
p+l,cr p-l,a 


-\2 * 
a u 


chi ^ 


+ 2kf 


p,a 


(19) 


Jp,a 

and forms the second derivative appearing in it according to the scheme 


r 

• 

8 u 



~ 8 u 

; 8 2 u* 



p,a+l 


p,a-l 


chi' 


p,a 


21 


( 20 ) 


from the first derivative 


8 u 


L^J 


already calculated in good approxi- 


P ,0 


mation according to equation (17) by jumping over. However, the case a = 1 
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Ik 


requires special consideration since 


du 


is not known at first. 


P,0 


But if one takes into consideration that according to equation (l6) 


u 




p,o 


dl 


one may put 


du 


Ibu 

_ 1 

d2ul 


^1 

P; 0 


P;1 


P,0 



5u 

+ 1 = 



Sq 

U -1 

pA 

d5 



du 


+1 


P,1 


- U Gp-l) 
2lS 


( 21 ) 


and hence calculate the value on the left from 
according to equation (17). 


ciu 




pA 


already known 


We now use the relation ( 19) not as a substitute for (15) in the 
sense that the entire calculation is to he made with ( 19 ); for it became 
clear - particularly near the separation point where the derivatives 

2ii become very large - that the convergence relations here can be easily 

blurred (unless an especially small k value was selected) ; the values 
u D obtained by iteration do not remain quite fixed; but creep on 
p ) a 

contlnously, although only by small amounts (compare also the remarks 
in section V ) . 

Rather we use equation (19) for making a correction in the values 
obtained after the last iteration in the manner described above. 


u l,0 

With the aid of the value 


<3u 




0 ,0 


(0 ^l); (already contained in the 


fourth column of table I) to which we add the value 


du 


dqj0;0 


just cal- 


culated according to equation (21) we determine (taking equations (20) 
and (14) into consideration) the values 


D 0 = u - 1,0 + ^ 


£^u 

* .if 

bu' 


du 

\ 

_ 

0 .a‘ - 1 '” 5 ( 


0,0+1 

Sq 

0,0-1) 


ill 1 I 1 ill ■ ill 
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We now assume, for Instance, that the values prescribed by 

table II were the final values even in the first procedure; we then 
Insert the values D CT in table II and calculate with the quantities ' 

A ct + B a appearing in them the values (corrected with respect to H,, 

|u 2^^ a = D a + 2(A 0 + B a ) 

these values, too, we note in the table. In the last column of this 
table we write the values 

= “ (tl) + Mi,, 

If the corrected values luoL deviate too much (that is, by more 

L Ji»cr 

than 1/4 to l/2 percent of U) from H a we calculate with the deri- 


vatives 


due 


hi 


Hi,, - H-1,1 


0,a 


2 k 


once more corrected values according to table III, p. 31 . The values 

?s] 1 ,, 

I = 


u, n then represent the final values for the profile at the point 
L j] 1,0 


For calculation of every step in 5 direction one must, therefore, 
calculate three to four of the calculation tables mentioned. The time 
expenditure may be estimated at approximately three to four hours per 
step. It should be stressed that all calculation operations are of 
purely schematic character and can therefore readily be performed by 
assistants. 

The values obtained are plotted on millimeter graph paper and the 
curve drawn through them. If slight scatter has resulted, after all, at 
one point or the other, one eliminates it with the aid of the drawing 
before starting on the next step. 

If the graph of the profile calculated just now shows that the curve, 
due to the increase in boundary-layer thickness, at the upper end no 
longer dies out gradually enough toward the asymptotic value U, one adds. 
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in calculating the following step, and q subdivision point in upward 
direction. 

The following condition should be mentioned which became evident in 
the practical calculation. If a step magnitude not sufficiently small 
is selected, two successive profiles may, due to accumulation of errors, 
show points where they are somewhat too close, or else somewhat too dis- 
tant from each other, compared to their actual course. In the calcula- 
tion this can be recognized by the fact that the third profile following 
these two profiles shows a behavior, at these points, compared to the 
second profile opposite to the behavior of the first compared to the 
second profile. For T) values at which the first two profiles were 
too close one notices a gap somewhat too wide between the last two and 
vice versa. If one does not want to repeat the calculation with smaller 
steps, one may, as was found practically to be useful, once omit the 
corrective calculation mentioned before for the profile to be calculated 
next, thereby eliminating the fluctuating of the profiles, and may then 
continue calculating in the normal manner. 


According to our calculation experiences one can recognize that the 
step magnitude k must be reduced in i direction by the fact that the 
two u values obtained in the corrective calculation which pertain to 


the same r 10 (thus in the example considered above the values 


and 


u 3 


L Jl,a 


l,a 


deviate from each other by considerably more than l/4 to 


l/2 .percent of the pertaining U value. 


If a new step magnitude in the £ direction, k]_, is selected kQ_< k 
^for instance, kq = there appears as a result, because of equation (l4), 
also a new step magnitude. 



in T) direction. If one wants to continue the calculation with the 
smaller steps kj_ for instance starting from I = I one needs as the 

initial values for further calculation the numbers 

u( !r->°^l) a( I p-k]_, ctZ]^) (ct = 1,2, . . .) 

The first named numbers may be read off directly on the profile curve 
for I = ! r already obtained. In order to obtain the latter, a double 
graphic interpolation must be made. One plots versus I the values 
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( 0 = 1^2 . . .) read off for the values of | = Ip^p = r^ 

from the curves of the previously calculated profiles. Generally it 
will he sufficient to do this for the values u^lp^oZi^ of three suc- 
cessive profiles, thus for p = r - 2, r - 1, r. From the curves drawn 
through them u = u(|,crZi) (a = 1,2 . . .) one may then read off the 
values u^ r -ki,crZi^ (a = 1,2 . . .). 

If a boundary layer is to be calculated up to the separation point’, 
it will in general be necessary to select, in the proximity of the sepa- 
ration point, rather small steps k. Since, however, due to the large 
increase in boundary- layer thickness, the profiles are here very elongated, 
one would obtain, because of the small step magnitude Z in q direc- 
tion, a very great number of subdivision points over the profile; this 
would of course increase the time expenditure for the calculation of a 
step. However, one may save a great deal of calculation expenditure by 
selecting, instead of equation (7)* for instance 


1 Su* _ 5 2 u* = 1 Su ji d ( u* , lAdu 1 dU(| ) 

2 si aq 2 u(£) St \J 0 ai \u(i) 2 ysT 2 


( 22 ) 


as the initial equation, and then performing the integration as before. 

If one again denotes the step magnitudes in £ direction by k, those 
in q direction by Z, one obtains instead of equation (l4) the relation 


k = 


Z 2 

k 


To the same Z as in the first considered case, therefore, there corre- 
sponds half the step magnitude in I -direction. 

In this manner the step reductions were carried out for the following 
examples of boundary-layer flow on the circular and elliptic cylinder. 

The convergence of the iterations now occurred no longer only on one aide 
toward the limit but alternately (except for the values assumed for small 
Tin* 


The numerical calculation showed further that a further step reduc- 
tion in | direction, still for the same Z, for instance with the aid 
of the initial relation 

i au* _ aV = _i_ au PV cm- ( + i\ au + 3 au(i) 

^ ai aq 2 u(5') aqJo ai \u(i) k 

was not advisable because the values assumed in the upper profile parts 
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on the right side are given as differences of two (approximately equal) 
large numbers and therefore scatter widely; by this the convergence rela- 
tions may he concealed. Thus, if one is forced to reduce the step mag- 
nitude k still further, one will do so in the manner described above 
with the aid of the relation (22) . It was found that one arrived in this 
manner, even for the extreme example of the circular cylinder, at a 
tolerable work expenditure even for the steps Immediately ahead of the 
separation point. 

The separation point I = l a (and therewith s = 
graphic interpolation, or extrapolation, of the values 
contained in the tables. 

The example of the Blasius flow at the flat plate shows very clearly 
•the high degree of accuracy attained with this method. Here the profile 
obtained by continuation could be compared with the exact profile. After 
calculation of six steps, the calculated values deviated so little from 
the exact ones that they could hardly be distinguished within the scope 
of drawing accuracy. The differences amount to less than l/2 percent 
referred to U. 

In order to enable following the mode of calculation in detail, we 
add the complete calculation of the first step in the continuation of a 
Blasius profile at the flat plate. 


is round oy 


8u 


IV. NUMERICAL EXAMPLES AND RESULTS 
1. Continuation of a Blasius Profile at the Flat Plate. 


The value s = 0 is to correspond to the leading edge of the plate. 
For the boundary- layer equation (3) which because of p'(s) = 0 Is 
simplified to 

ha + - = 

8s ^ 8 t| 3q 2 

together with the continuity equation 


+ 8V = 0 
8s 
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then, exists according to Prandtl-Blasius, as is -well known, a solution 
of the form 


v s = J Up’(l) with 


for which applies 


v s — VO for 



tj— VO 


and 



s — VO and all T) 


v s - 


>u 


for 


and all s 


^ 0 


The function cp = 9(0 satisfies the ordinary differential equa- 
tion of the third order 


9 9' 1 = - 9* ' ' 


and the boundary conditions 

9 ( 0 ) = 0 9 ' ( 0 ) =0 lim 9’(0 = 2 . 

£ — > 00 


The values of 9 ’ ( C ) are to he taken from the following table: 


i 

\ <p'(0 

0 

0 

0.1 

0.0664 

0.2 

0.1328 

0.3 

0.1989 

0.4 

0.2647 

0.5 

0.3298 

0.6 

0.3938 

0.7 

0.4563 

0.8 

0.5168 

0.9 

0.5748 

1.0 

0.6298 


i 

\ <p'(0 

1.1 

0.6813 

1.2 

0.7290 

1.3 

0.7725 

1.4 

0.8115 

1.5 

0.8460 

1.6 

0.8761 

1.7 

0.9018 

1.8 

0.9233 

1.9 

0.9411 

2.0 

0.9555 


5 

| 9'(0 

2.1 

0.9670 

2.2 

0.9759 

2.3 

O.9827 

2.4 

O.9878 

2.5 

0.9915 

2.6 

0.9942 

2.7 

0.9962 

2.8 

0.9975 

2.9 

0.9984 

3.0 

0.9990 


We choose U = 1 so that we may put 


| = s 
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and start our continuation procedure at s =. 1. As step magnitude in 
V direction we take l = 0.6 so that k becomes equal to 0.18. The 
initial profile then may be taken directly from the above table, whereas 
the profile one step farther back, thus the profile at s = 0 . 82 , is to 
be obtained from this table by graphic interpolation. The values are 
contained in the following table: 


7 

u( 0 . 82 ,T]) 

0 

0 

0.6 

0.225 

1.2 

0.430 

1.8 

0.626 

2 .4 

0.782 

3.0 

0.894 

3.6 

0.955 

4.2 

0.982 

4.8 

0.995 

5.4 

1 


7 

71(1,7) 

0 

0 

0.6 

0.1989 

1.2 

0.3938 

1.8 

0.5748 

2.4 

O.7290 

3.0 

o. 846 o 

3.6 

0.9233 

4.2 

0.9670 

4.8 

0.9878 

5.4 

0.9962 

6.0 

0.9990 

6. 6 

1 


Six steps (that is, up to s = 2.08) were calculated by the method 
described. The calculation of the first step is contained completely in 
the table added at the end of the report. The results are represented 
in figure 2. 


2. Circular Cylinder According to Hiemenz. 

Hiemenz 9 measured the pressure distribution on a circular cylinder 
of diameter 2r = 9.75 centimeters i m mersed in water and approached by 
the flow at a velocity of 19.2 centimeters per second. % 

In order to make the quantities appearing in the basic equations 
dimensionless, one introduces the reference length 1 = 1 centimeter 
and the reference velocity V Q = 7*151 centimeters per second which corre- 
sponds for v = 0.01 centimeter^ per second to a Reynolds number 

7V 

K = -V* “ 715.1 


%. Hiemenz: Die Grenzschicht an einem in den gleichformigen Fluss- 

igkeitsstrom eingetauchten geraden Kreiszylinder (The boundary layer on 
a rectilinear circular cylinder immersed in the uniform fluid flow). 
Dissertation Gottingen, 1911, published in Dingier 's polytechn. J. Vol. 326, 
1911, pp. 321-342. : 
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then, the velocity distribution measured for 0 = s ='1 } that is up to 
the separation point, observed shortly before s = 7 (corresponding to 
an angle a of 80 ° to 82° from the forward stagnation point) may be 
represented satisfactorily by the polynomial 


U(s) = s - 0.006289 s3 _ 0.000046 s5 

On the basis of the previous indication, the solution of Blasius- 
Hiemenz could be used up to the value s = 4. 5 (a ~ 55°) so that our 
calculation starts at s = 4.5 ( a s does Gortler ' s-^ ) . The value l = 0.4 
and thus k = 0.08 were selected as step magnitudes for the first steps. 
The representation of 

1 = /J 5 iW dt and 0 ' u(s) 


against s may be seen from figure 3. 


The initial profiles at 4 = - 0.08 and 
tables, are compiled, together with the values 
from the same tables, in the following table: 


n 

u(-0,08,tj) 

u(0,T| ) 

8v s 

8s 

s=4.5 

0 

0 

0 

0 


0.4 

1.282 

1.289 

0.008 

0.8 

2.229 

2.2268 

0.104 

1.2 

2.871 

2.948 

0.249 

1.6 

3.269 

3.384 

0.367 

2.0 

3.488 

3.628 

0.448 

2.4 

3.'596 

3.749 

0.497 

2.8 

3.649 

3.813 

0.529 

3.2 

3.667 

3.834 

0.538 

3.6 

3.672 

3.839 

0.540 

00 

3.674 

3.842 

1 



4 =0, taken from Howarth's 
8u 


of 


3=4.5 


resulting 


When the latter values are used, the calculation of the first step 
requires only one worksheet of the type described before. With the step 
magnitudes indicated) first four steps (up to 4 = 0 . 32 ) were calculated. 


10 


See footnote 3* 
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The profiles obtained are represented together with the initial profiles 
in figure 4 (partly displaced with respect to each other). 


With the initial relations (22) as a basis, five further steps 
(up to I = 0.52) were calculated for the same Z = 0.4 and the required 
k = O.O^. Like wise with the use of equation (22), one step (| = 0.5^) 

with Z = \fo.08 = 0.283 and k = 0.02 and finally two more steps 
with Z = 0.2, k = 0.01 (up to I = 0.5 6) were calculated. The pro- 
files are also represented in figure 5* 


Sul 


Pp> qI ) 


= 6.87 


the separation point was found 
( compare figure 6) . Thus 


By plotting of the values 
to be | sep = 0.5697, that is s sep 

all together twelve steps were to be calculated. Figure 7 shows the 
curves u = u(| , const). Their steep decline in the neighborhood of the 
separation point is remarkable. 


Figure 8 shows a comparison of a few of the profiles obtained by us 
( -S) with those of Blasius-Hiemenz ( — — B-H), Pohlhausen ( — - — P), and 

Gortler ( G) which were obtained for the same pressure distribution 1 !. 

The comparison shows, first of all, that the Blasius-Hiemenz solution 
becomes insufficient in the neighborhood of the separation point; the 
reason obviously lies in the fact that the series developments used con- 
verge for large s only slowly, if at all, and that, therefore, with 
merely the first three terms the actual course is not satisfactorily 
represented there. 


Our values agree best and most systematically with those obtained 
by Gortler. The differences are increasingly noticeable toward the sepa- 
ration point. The deviations from the values obtained by Pohlhausen, 
considered as a whole, remain for this example within tolerable limits 
although a systematic variation of the differences cannot be determined. 
It is remarkable that the differences assume higher values precisely in 
the proximity of the velocity maximum (| ~ O.36, s ~ 6, a ~ 71°) (com- 
pare the curve for | a 0.32 represented in figure 8) while again sub- 
siding to some extent toward the separation point. 

The separation point was found according to Gortler in good agree- 
ment with our value s se p = 6.8, according to Hiemenz at s 8ep = 6.98, 
and according to Pohlhausen at s se p = 6.9^. An approximately correct 

position of the separation point is, therefore, by itself not yet deci- 
sive for the usefulness of a method. 

HCompare K. Hiemenz, paper quoted in footnote 9, H. Gortler, paper 
quoted in footnote 3, and K. Pohlhausen, "Zur naherungsweisen Integration 
der Differentialgleichung der laminaren Grenz schicht" (On the approximate 
integration of the differential equation of the laminar boundary layer), 
Z.A.M.M. Bd. 1, 1921, pp. 252-268. 
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3 . Elliptic Cylinder of the Aspect Ratio 1 : 4 . 

As a further example, we calculated the boundary layer for an ellip- 
tic cylinder of the aspect ratio 1 : 4 , taking as a basis the pressure 
distribution resulting from the potential theory. 


I = |g. and V 0 = 4.3 U Q 

were chosen as reference quantities for the introduction of dimensionless 
quantities, with l Q being half the circumference of the ellipse and 
U Q the free stream velocity. The (dimensionless) velocity at the edge 
of the boundary layer could be taken directly from a table by Schlichting 
and Ulrich.-^ It is represented in figure 9 together with the function 


£ = i (s) 


/ ft 

Jo. 2 uTtT 


In the interval 0 % s = 0.2 
satisfactorily by the polynomial 


it was possible to represent U = U(s) 


U(s) = s - 5.116 s^ 


The initial profile, however, was chosen at s = O.163 (£ = - 0 . 25 ) for 
the reasons mentioned before. For the first six steps l = 0.5 and 
k = 0.125 were taken as step magnitudes. The two initial profiles are 
represented in the following table, together with the values 
(Bvs^ 



0.163 


71 

u( -0.375 ,ti) 

u(-0.25,T)) 


8v s 

_8s 

s = 0.163 

0 

0 

0 

0 


0.5 

0.0573 

O.O589 

0.0949 

1.0 

0.0951 

0.0990 

0.2672 

1.5 

0.1166 

0.1226 

0.4235 

2.0 

0.1268 

0.1340 

0.5197 

2.5 

0.1310 

0.1389 

0.5703 

3.0 

0.1323 

0.1403 

O.5858 ■ 


0.1327 

0.1408 

1 



12 h. Schlichting und A. Ulrich, "Zur Berechnung des Umschlages lam 
turbulent" (On the calculation of the transition from laminar to turbu- 
lent) Bericht S 10 der Lilienthal-Gesellschaft ( 1940 ), pp. 75-135- 
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From | = 0.5 onward the steps could, first, he increased. The 
following further steps were calculated: Seven steps with 

2 = \jo. 5 = 0.7071 and. k = 0.25 (up to I = 2.25), three steps with 

2 = 1 , k = 0.5 (up to i = 3 - 75 ), four steps with 2 = l/i = 1.4l4 
and k = 1 (up to | = 7 . 75 ), and nine steps with 2=2 and k = 2 
(up to | = 25.75). 

With the aid of the starting equation (22) another step reduction 
was made. The further steps were: Three steps with 2=2 and k = 1 

(up to I = 28 . 75 ) and finally one step with 2 = \/2 = 1.4l4 and k = 0 . 
(up to I = 29.25). A complete calculation was thus made of 33 steps 
altogether. 


It became clear that selection of larger steps is not advisable, 
particularly at the point where the curve U = U(s) turns from its 
steep ascent to the flatter course (compare figure 12 ). 


On the "high plateau" of velocity distribution itself one could 
have chosen steps somewhat larger but they would have had to be reduced 
again when approaching the separation point. A large number of the pro- 
files we calculated can be seen in figure 10. The separation point was 


determined from the variation of 
figure 11 ). 



as s se p = 8.475 (compare 


Since Schlichting and Ulrich completely calculated^ the same 
example once according to the ordinary Pohlhausen method (P^-method) , 
and then according to a Pohlhausen method modified by taking a polynomial 
of the sixth degree as a basis (Pg-method), the comparison could be made 
for a number of profiles. The results are compiled in figures 13 and l4. 


As far as the pressure minimum the deviations between our curves 
and the P 4 - and Pg -curves are not too large. However, larger deviations 
appear in the proximity of the separation point. There the profiles of 
the original Pohlhausen method agree with ours better than the profiles 
of the Pg-method, especially for small q values. The resulting separa- 
tion point was, according to the P 4 -method, at s ge p = 8 . 38 , according 
to the Pg-method, at s ge p = 8 . 26 ; thus these values (especially that 
of the P^ -method) do not deviate too widely from our value. 


■^Quoted in footnote 12 . 
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V. REMARKS REGARDING THE CONVERGENCE OF THE ITERATION PROCESS 


Regarding the conditions of convergence of the iteration process 
described in section III, ve can prove the following theorem which will 
probably be sufficient for the requirements of practical calculation. 

If there applies for the profile at the point | ='| p for all =j t}CT 0 


° " U(lp) 

^p ,1 Il <1 + 

U(V) u(5p) 

the sequence of the velocity values obtained by the iteration process 


and 


chi 


[ U nlp+1^CT - a o) 


converges with increasing n. 

Thus one is always able to predict, when calculating a new step, 
whether the iteration process will converge. The presuppositions of 
the theorem are satisfied with certainty when for all T] ^ qa 0 

0= [u^ ^ < U(lp.) and 

apply as is the case for instance for the profiles before the pressure 
minimum; for then 




^2 


^ 0 




3u 

< 

= Ti 

5u 

ill- 

^ 0 ^. ' 





= ,n < u Up) 


is selected so that 




du 




(tj = T)cr 0 ) 


is valid if 
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The presupposition's are satisfied also for the S-shaped profiles 
beyond the pressure minimum when does not become too large. For 

the examples we calculated, boundary-layer flow for circular and elliptic 
cylinder, the presuppositions are satisfied, up to the separation point. 
Proof of the theorem: If one takes into consideration that according 

to the procedure of section III one has to put 


* 

u n+1 


u *p,a+l + u *p,a-l 


P+1, a 


U(ip) 



( 

Su 

l \ 

-^1- 

< 


p+l/r “ u P-l/r 


T=1 


2k 


1 

2 



k 

U(l p ) 




p+l,g ~ 


2k 


there follows with 


r - 

r 

r 





tin 

= Max - 

** * T=l,2,...,a 

u n 

P+1,T 

u n-l 

P+ 1 ,T 



with the presuppositions taken into consideration, obviously 


uufl 


P+1,1 


u u 


P+1,1 


< 1 


2U(ip) x 


chu 


^pjTil 


Wp+1,< 


2U(^p) 


( TT ( | p) - Hpn,a 


< „ 
= a 


Hcti,. 
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with 0 < a < 1 and a beind independent of n. Furthermore, one can 
see for oneself that a (0 < a < 1, independent of n) can be chosen 
so that simultaneously the estimations 

[ Uu+ 1 p+ 1,0 " Uu P+l^ = a d u p+1, a 


• a d.u p+i^a 

exist so that 

dn+1 n = a d n 
u+ p+1, a 11 p+1, 0 

must be true as well. 

Hence there exists the limit 

^1,0 = lim ujx 

T)— +0O |_ J 

U1 0.1,0 + ( U2 0.1,0 ■ U 1 p»l,o )* + 

([““Jo.1,0 - [ u n-l]p.l,o) + •" (o j 0 °> 

since the series at right may be majorized by a convergent geometrical 
series. 

One recognizes further that - if equations (l8) or (19) are used 
instead of equation ( 15 ) - if convergence of the iteration process 
can be proved under the assumptions that for | = 4 P and all r\ = 

0 S M ep> , < »<tp> 



I"™- 1 ] 0 . 1,0 


u P+1, a 
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and 


du 


< !>]* „ 


since one then finds the estimations 

n, U(!p) - [u] 


u u+l 


P+1,1 




•Jp+1,1 


< 


U(6n) 


*0^1 


UUo) 


hi 


p+l, a 


u u+l 


PH-1, cr 


% 


P+1,0 


' A 5(lp) ' r " J e p ^ c 


U(lc) 


U(£ 0 ) 


[ d Jp+l, a 


thus 


hi+l 


P+1,0 


< 

= a 


x u 


P+1, a 


again with 0 < a < 1 and a independent of n. These presuppositions 
are satisfied for instance for the profiles before the pressure minimum. 
However, in the proximity of the wall, if the profiles there show an 
approximately rectilinear course, the convergence will take place only 
very slowly. 

Beyond the pressure minimum one can, therefore, not arrive at a 
general statement on the convergence. As mentioned before, our calcu- 
lations in the proximity of the separation point showed that the case 
of divergence may actually occur. Therewith the procedure we selected,, 
using the relation ( 19 ) merely for the correction calculation, proves 
to be perfectly reasonable also from the general point of view now 
considered. 

Correspondingly, one recognizes that the iteration process performed 
with the aid of relation (22) certainly is convergent at the point 
I = |p for 0 ^ 0 O if there for all tj s t}0 o 


0 - 
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If one takes generally the initial relation 


1 Su* S^u* 

m S| St}2 


_i_ r Su dT1 
u(i) on J o si 



Su 

SI 


+ 



convergence at the point | = | p for a = a Q would he assured under 

the following sufficient conditions: For all t] = T]cr 0 there shall he 
valid 


and 




Ill I I I IHII HUB 
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Thus, if one wants on the right side a positive limit also for 

U(6 p ) 


m must he not greater than 3. * 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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TABLE I. 



TABLE II. 


b'l,. 

- 1- Bq + C 0 

Mi,* 

■ M.,. * ' <E ‘> 

% 

[*% 

= D(j + 2(A a + Bg) 

r R, 

■ * '<>»> 

-U(Si) 

* 

0 

* 

* 

-U(5 X ) 

* 

0 

* 

* 

* 

* 

* 

* 


TABLE III, 
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THE FIRST STEP IN THE CONTINUATION C 

JF 


THE BLASIUS PROFILE 



I 




A a 


Ca 

* * 



= -<%e<r 

a<j + Ba . u °^; u v" 

® + B(j + C ( 


-0.1450 
-0.'346l 

- 0.2011 
-0.4855 

-0.2844 

- 0.5788 

2.4 -0.2944 

- 0.5611 

3.0 -0.2667 

-0.4428 

3.6 -0.1761 

-0.2594 


o. 05027 
0.04068 


-0.14420 
-0.10912 
-0.07654 
8 


-0.02772 

-O.OI381 

-0.00594 

- 0.00220 

-0.00068 


-0.0209 -0.0235 
■0.0219 -0.0302 
-0.0218 -0.0365 

-0.0144 -0.0334 

-0.0259 
-0.0024 -O.OI63 

-0.0005 -0.0087 

-0.0001 -0.0040 


.8031 

.6132 

.4386 

.2896 

-0.1739 

-0.0935 


- 1.0000 

-0.8266 


-0.0066 

-0.0019 

-0.0005 


-0.0532 

-0.0224 

- 0.0081 

-0.0019 

-0.0005 
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Figure 2.- Continuation of a Blasius profile at the flat plate. 



4.4 



Figure 3.- The functions | = g(s), U = U(s) for the circular cylinder. 










Figure 7.- The curves u = u(| const.) for the circular cylinder. 
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Figure 8.- Comparison of the profiles for the circular cylinder. For 
I = 0.16 the Blasius-Hiemenz profile coincides partly with Gortler’s 
and partly with our profile. 
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Figure 9(b).- The functions I = i(s) and U = U(s) for the elliptic 

cylinder. 





Figure 10.- Velocity profiles for the elliptic cylinder. 
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